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We study the dynamics of phantom model and give the conditions for potentials to admit tracking 
attractor, de Sitter attractor and big rip attractor. Especially, we show that phantom models with 
exponential and inverse power law potentials do not admit the tracking attractor solution. The 
"tachyonic" stability of the system at/near the attractors as well as the quantum stabiltiy have 
been studied. 

PACS numbers: 98.80.-k, 95.35. +d 
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I. INTRODUCTION 

Astronomical observations on the CMB anisotropy 
P, 0, 0, the relation between red-shift and luminosity 
distance of Supernovae 0, 0, as well as the galaxy dis- 
tribution, clumping and spreading depicted that our 
Universe is spatially flat, with about two-thirds of en- 
ergy density resulted from dark energy that has an equa- 
tion of state w < —1/3 and accelerates the expansion 
of the Universe. The nature of this substance is quite 
unusual and there is no justification for assuming it re- 
semble known forms of matter or energy. Candidates for 
dark energy have been widely studied and focus on a dy- 
namically evolving scalar field (quintessence |3.l9l lTol[lll |. 
with w > —1 and phantom with w < —1) and cos- 
mological constant. Present observation data constrain 
the the range of the equation of state of dark energy as 
— 1.38 < w < —0.82 [13, which indicates the possibihty 
of dark energy with w < —1, debuted as Phantom or 
super quintessence . The realization of tx; < — 1 could 
not be achieved by scalar field with positive kinetic en- 
ergy and thus the negative kinetic energy is introduced 
although it violates some well known energy conditions 
. Another important consequence of Phantom is the 
Big rip |1^] or Big smash phase, in which the scale 
factor of the Universe goes to infinity at a finite cosmo- 
logical time. The cosmological implications of Phantom 
have been widely studied 17, 18, 19, 20, 2i, ^ ^ .24, 
mmmiiillallilllllil ,thc Phantom model with 
Born-Infeld type Lagrangian has been proposed '33| and 
its generalization to brane world has been done in Ref. 

One of the most important issues for dark energy mod- 
els is the fine tuning problem and a good model should 
limit the fine tuning as much as possible. The dynamical 
attractor of the cosmological system has been employed 
to make the late time behaviors of the model insensitive 
to the initial condition of the field and thus alleviates 
the fine tuning problem. In quintessence models, the dy- 
namical system has tracking attractor that makes the 



quintessence evolve by tracking the equation of state of 
the background cosmolo gica l fluid so as to alleviating the 
fine tuning Droblem|3(l l37|| . While in phantom models 
with canonical lagrangian, the kinetic energy term be- 
comes negative, which make the phase space of the cos- 
mological system behaves rather differently from those of 
the quintessence. The major difference is that the track- 
ing attractor does not exist in the phantom system with 
exponential and inverse power law potentials, which ad- 
mit tracking attractor solution in quintessence models. 
Accordingly, We give the condition for the phantom sys- 
tem to admit tracking attractor. On the other hand, 
there are also two late time attractors in the phantom 
system corresponding to the big rip phasejs^ and de Sit- 
ter phasepOj. One cannot say, as a priori, whether our 
Universe will end with the big rip or with the de Sitter 
phase, therefore it will be interesting to study dynamical 
evolution of the Phantom models in general potentials 
and the conditions for big rip and de Sitter phase respec- 
tively. 



II. BIG RIP AND DE SITTER LATE TIME 
ATTRACTORS 

Since current observations favor flat Universe, we will 
work in the spatially flat Robertson- Walker metric. The 
corresponding equations of motion and Einstein equa- 
tions could be written as. 



H = -{p-y+Pj - r) 

P7 = -3i?(p7 +P7) 
4> + 3H^-V'{(t)) =0 

^2 



(1) 
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where — 8ttG, is the density of fluid with a bary- 
otropic equation of state — (7 — l)p-y, where < 7 < 2 
is a constant that relates to the equation of state by 
ui = 7 — 1; The over dot represents derivative with re- 
spect to t, the prime denotes derivative with respect to 
(f). Pel, = -^0^ + V{(l)) and p^ = - V{(j)) are the 
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energy density and pressure of the (j) field respectively, 
and H is Hubble parameter. 

However, for an arbitrary potential, and considering 
the presence of other energy density, one cannot find 
analytically solvable models. Thus, we need to ana- 
lyze the models via phase space analysis. Similar as in 
Ref. [3^ I37I Ib^ , we introduce the following dimension- 
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less variables 



V6H 



y = 



V3H 



' kV(0)' 



r — ^ y)]^^!'^'' and N = log a. Then, the equation sys- 
tem (^1 could be reexpressed as the following system of 
equations: 



Also, we have a constraint equation 
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The equation of state for the scalar fields could be ex- 
pressed in term of the new variables as 
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X — y 



(4) 



Different from the case in exponential potential |39l |. 
the parameters A and F here are variables dependent on 
(j). Thus, strictly speaking, the above system is not an 
autonomous system. So, just as was done in Refs. 
l37j . if we want to discuss the phase plane, we need certain 
constraints on the potential, or equivalently the condition 
under which the potential has the property we want, in 
order that we can get some explicit results. 

In the following, we will read out the big rip attractor 
as well as de Sitter attractor from Eqs.(|2) and specify the 
corresponding conditions for the potential. The big rip 
attractor is the dynamical attractor that corresponds to 
the Phantom domination — 1 and an equation of state 
w < —1. While, the de Sitter attractor also corresponds 
to Phantom domination but with w = —1. Looking at 
Eqs.Q, one can observe that physically meaningful criti- 
cal points (xc, yc,Xc) of the system are (Note that here we 
are only interested in the expansionist solutions {y > 0)) 

(i) . (-^, V'l + ^, Ac) with F = 1 and Ac ^ 

(ii) . (0,1,0) for all F. 

(iii) . (0,0, Ac) for all F. 

where Ac could be fixed by F = 1. To gain some in- 
sight into the property of the critical points, we write 
the variables near the critical points (a;c,yc,Ac) in the 
form X = Xc + u, y = yc + v, and A = Ac -l- <^ with 



u, V, S the perturbations of the variables near the crit- 
ical points. Substitute the expression into the system 
of equations (|2Jl, one can obtain the equations for the 
perturbations. The coefficients of the perturbation equa- 
tions form a 3 X 3 matrix M whose eigenvalues determine 
the type and stability of the critical points. 

(i) Big rip attractor: in the above case (i), when F ~ 1, 
from Eq.(|2Jl, we know that A is almost constant. So, sim- 
ilar as was done in Refs. [s^ [s^l, we consider only the 
first two equations of Eqs.Q. Thus, the corresponding 
matrix M should be two dimensional and the correspond- 

ing eigenvalues are (—3 — ^, —37 — Ac), which indicate 
that this critical point is a dynamical attractor of the 
system. It is not difficult to evaluate that at the critical 
point, the cosmic energy density parameter fi^ = 1 and 

the equation of state w = —1 f. So, this attractor 

corresponds to an equation of state w < — 1 and will lead 
to the catastrophic big rip. It is worth noting that for an 
arbitrary potential, Ac is not a constant, which could be 
determined by solving F = 1 for a given potential. While 
in the exponential potential, Ac is denoted as A and is a 
constant parameter of the potential. For a general po- 
tential V{(j>), the condition for the existence of the above 
mentioned big rip attractor is F = l(for practical pur- 
pose, this condition could be F ~ 1, as was done in Refs. 
IMIS) and X^O. 

(ii) de Sitter attractor: in the above case (ii), the corre- 
sponding eigenvalues of matrix M are (—3, 0, —37), which 
indicate that the critical point is a dynamical attractor 
corresponding to equation of state w = — 1 and cosmic 
energy density parameter fi^ = 1, which is a de Sitter 
attractor The condition of such a de Sitter attractor 
is A = 0. According to the previous definition for A, this 
condition is equivalent to that the attractor appears at 
the extremum of the potential. We must point out that 
when keeping the perturbation to the linear order, the 
right hand side of the third equation in Eqs.Q will van- 
ish. Thus the condition A = is only a less restrictive 
one, from which we cannot determine if the extremum 
is maximum or minimum. To extract more information, 
one needs to either keep higher order perturbation or 
choose different dimensionless variables as. was done in 
Ref. li^l , in which we choose x = j^, y = and reduce 

Eqs.lpJ to a different form. The analysis indicates that 
the extremum must be a maximum for the existence of de 
Sitter attractor for Phantom models. For details, please 
refer to Ref. 

In the above case (iii), the critical point corresponds to 
— and the corresponding definition of becomes 
meaningless. Moreover, the corresponding eigenvalues of 
the M matrix is (0, 3(7 — 2)/2, 37/2), which indicates the 
critical point is not a dynamical attractor. So, we won't 
address it into details. 

Next, let's consider two specific models that admit big 
rip attractor and de Sitter attractor respectively. To do 
so, we must specify the potential. For the big rip at- 
tractor, we have the requirement that A 7^ and F = 1 
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FIG. 1: The evolution of equation of state w vs. A'^. We 
choose dimensionless parameter — — 1.2 and the j = 1. 



FIG. 3: The evolution of equation of state w vs. TV, in which 
for simplicity, we choose Vb = 00 = k = ao = 7 = 1 
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FIG. 2: Phase graph of the model for different initial x and y. 
We choose dimensionless parameter — = 1.2 and the 7 = 1. 
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FIG. 4: The phase graph of the system with x = (j)/4>o and 
y — 4>/cj>o, which are so defined for convenience of evaluation. 
The condition for evaluation is Vo = (j>o = k = ao = ■y = 1. 



(or r ~ 1 as approximation for practical purpose). It 
won't be difficult to find that the solution for the exact 
r = 1 would be that = Cie^^'^ with Ci and C2 as 

arbitrary constants. Conventional study for quintessence 
models focus on the case that C2 < 0. Here, C2 could 
be both positive and negative constant. For definite, we 
choose V{(j)) = Voe°"^ with a a positive constant, and 
thus A — In Fig.l and Fig. 2, we plot the numerical 
results. 

For the de Sitter attractor, we require A = 0, which is 
equivalent to that the critical point is the extremum of 
the potential. We choose the potential 



V{cj,)^Vo-ao{^)' (5) 



as a toy model, in which the de Sitter attractor corre- 
sponds to the maximum of the potential Eq.(|SJ|. The 
corresponding numerical results listed in Fig. 3 and Fig. 4 



III. TRACKING ATTRACTOR 

In this section, we will show the tracking attractor be- 
havior of the dynamical system that are not manifested 
in the above analysis. We consider the situation that A 
is very large and F is not 1 but nearly constant. It will 
be convenient for the discussion if we make the following 
transformation e ~ j, x — eX and y = ey|37j. Then the 
dynamical system could be rewritten in terms of the new 
variables X and Y as: 

^ ^_V6(r-i)x2 + ^7X-3x-y|r2 (6) 

— =-V6{r-l)XY- \[^XY+-jY 
dN ^ ' V 2 2 ' 

where we have omitted the terms contain e. The 
critical points of the above autonomous system Eq.|(Bl) 
corresponding to the tracking solution is {Xc,Yc) — 

\J%le{le ~ 2)), whcrc 7e = & defined 

effective barotropic index of the scalar field. The cor- 
responding cosmic energy density parameter for phan- 
tom is fi^ = and the equation of state for phan- 
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torn is W0 = with = 7e — 1. Clearly, to make 
il^ physically meaningful, we need F < 1/2 and there- 
fore 7e < 0. Since A is large, then fi^ is small and 
as a consequence, the background dominates. Also, at 
this point, the energy density parameter of the phantom 
field will track the effective barotropic index of the scalar 
field 7e, by which the name of tracking solution is jus- 
tified. It is not difficult to observe that although the 
equation of state of the background cosmological fluid is 
positive (0 for matter and 1/3 for radiation), the effective 
equation of state could be less than -1, which is required 
for the the cosmic density parameter for phantom to 
make sense. This tracking property is very different from 
that of quintessence models in which the corresponding 

critical point is {Xc,Yc) = ( y^Te, y^f 7e(2 - 7e)) with 
the 7e defined the same as that in the above phantom 
case. But the corresponding cosmic density parameter 
is Ctquint = ^ which is just the opposite of the ex- 
pression for phantom. This sign difference lead to dif- 
ferent requirements for F, that is, for quintessence, F 
must be greater than 1/2 while for phantom F must be 
less than 1/2 so as to make their respective cosmic den- 
sity parameters physically meaningful. This constraint 
on the potentials for phantom excludes the exponen- 
tial potential(F((^) = Voe°"^ with F = l)and inverse 
power law potcntial(y((/)) = ^/^" with F = 1 + i) that 
have been used in quintessence models. A simple po- 
tential that satisfies the above condition takes the form 
V{0) = A(f)^+T^ with n> 1. In this case, F = j^. 

Next, let's show the stability of the critical point. To 
do so, we use the similar approach used in previous sec- 
tion and find out the eigenvalues of the perturbation 
equation of the dynamical system near the critical point. 
We write theX = Xc + U and Y = Yc + V, with U, V 
the perturbations. Then the equation system could be 
written as, up to the first order of the perturbations. 



dU 

dN 

dV 

dN 



= [-3 + (9/2 - 3F)7e]C/ - 3v/7e(7e - 2)V (7) 



-3/2(2F-l)y^ 
The eigenvalues are 



2)U 



Ai =--|2 + (2F-3)7. 



(8) 



+ V4 + (4 - 24F)7e + (1 + 2F)272 



2 + (2F - 3)7e 



-v/4 + (4 - 24r)7e + (1 + 2r)272 



-1.5F 



-2 



V -2.5 



-3.5 



0.1 0.2 0.3 0.4 0.5 

Gamma 



FIG. 5: The second eigenvalue A2 vs. F at the matter domi- 
nated epoch. 
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FIG. 6: The second eigenvalue A2 vs. F at the radiation 
dominated epoch. 



second one, the term in the square root become negative 
if r < ^~^^ when 7 = 1 (matter dominated epoch) and 

P < when 7 = 4/3 (radiation dominated epoch). 

In these two cases, the critical point is stable. While 
for the range of < F < 1/2 and < F < 

1/2 corresponding to the matter and radiation dominated 
epoch respectively, we plot the second eigenvalue vs. F 
in Fig. 5 and Fig. 6. From the plots, one can see that 
both cases yield negative eigenvalues and therefore the 
critical point is stable for the range F < 1/2. That is, the 
tracking solution is a dynamical attractor of the system. 



IV. TACHYONIC STABILITY NEAR THE 
ATTRACTORS 

In this section, we will study the "tachyonic" instabil- 
ity^ of the dynamical system near the dynamical attrac- 



According to our previous discussion, we know that F < 
1/2 is required for the solution physically meaningful. 

Closely examine the two eigenvalues, one can find that 
the first one always has a negative real part. While for the 



^ Wc call the system is tachyonically instable if the effective mass 
of the perturbation becomes imaginary 
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tors. The reason for us to discuss the stabihty only near 
the attractors is that we consider that our universe is just 
at or near one of the attractors. The stabihty of the sys- 
tem at the attractors does not necessarily guarantee the 
stability during its evolution, but if we took a anthropic 
philosophy and the above assumption, we can say that we 
are just living in a universe where the phantom evolves 
safely to its stable attractors. In synchronous gauge, the 
perturbed equation of motion reads [13. IT^ 



+ iHScPk + (k^ - V"{<P)) = (9) 

where the h is the trace of hij (the metric perturbation) 
and S(f>k is a fourier mode of the phantom field perturba- 
tion. The prime denotes the derivative with respect to 
(j) and the dot denotes the derivative with respect to t. 
The effective mass for the perturbation is [k'^ — V" {4>)]^^^ . 
For the potentials with negative V" {(/)), it will be always 
stable. While for potentials with positive y',the stabil- 
ity requires that the wave number of the perturbation 
should be fc < kcrit — \/V"{(p). In terms of the dimen- 
sionless variables we introduced previously, the critical 
wave number could be written as 



V. QUANTUM STABILITY 

Besides the tachyonic instability we have analyzed in 
section IV, there is also another instability quantum in- 
stability, that plagued phantom models Since dark 
energy is assumed to be weakly interacted with other 
matter except via gravity, it would be natural to con- 
sider only the phantom decay relevant to gravitons, i.e. 
a phantom decays into other phantoms and gravitons. 
Clearly, the simplest process will be one phantom decays 
into two other phantoms and a graviton as 

0,^/1-1-01+02 (11) 

Follow the discussions in Ref. we firstly consider 

the decay in the tree level with potential coupling and 
calculate the corresponding Xeff for the potentials in 
this paper, i.e. exponential potential 1^(0) = Voe°"^, 
quadratic potential Eq .(O and an analytically solvable 
potential(its analysis will be presented in section VI). 

V{(f>) = Vn{sm[V3K{(j> ~ 0o)] -f csc[V3k{<P ~ 0o)]} (12) 

It is not difficult to evaluate that the corresponding 
Xeff are 



x/l^"(0) = V^Hy\X\Vf 



(10) 



At the tracking attractor, we can easily obtain that 
Uc = \ ^"''^^^"a"^^ and the corresponding critical wave 



number will be k^rit — 3iJ 



r7e(7e-2) 



If we choose the 



Hubble parameter H as the present value and 7e as well 
as r are of order 1, then the stability will appear if the 
wave length of the perturbation is greater than about 
10^* cm, which is even greater than the radius of our ob- 
servable Universe. That is to say, the tracking attractor 
solution corresponds to a state that is "tachyonically" 
stable. At the big rip attractor, the critical wave number 



should be kcrit — \Xc\H^'i+ with a dimension- 
less variable that can be specified by solving F = 1. If we 
also choose the Hubble parameter H as the present value 
and Ac is of order 1, then the situation will be similar 
to the above case for tracking attractor and the model is 
stable for practical purpose at the big rip attractor. At 
the de Sitter attractor, we have V"{(t)crit) < so that 
the point corresponds to the maximum of the potential 
[40l |. Therefore, the models will always be stable at the 
de Sitter attractor. From above discussion, we know that 
the de Sitter phase, big rip phase and tracking phase all 
could not only be dynamical attractors under properly 
constrained potentials but also stable state in terms of 
the tachyonic stability. However, as we have mentioned 
at the beginning of this section, the stability at the at- 
tractors does not guarantee its evolution. But this could 
be circumvented by adopting an anthropic view. 



V// 



" gm'"^ ' exponential potential 
, for quadratic potential 

for potential Eq.lO 



(13) 



Similar as in Ref.0|, 0o is of the order of Mp. Note 
that the parameters in the potentials are of the order 
a K ^ 1/Mp, the Xeff for exponential potential as 
well as the solvable potential (|12|) are both of the order 
Vo/A/p ^ 10^12", the same as that in Ref. 14] and thus 
need the same momentum cutoff. But for the quadratic 
potential, Xeff = and therefore the corresponding 
model need no momentum cutoff and is stable in the 
tree level potential coupling. As we have clarified at the 
beginning of this section, the decay described above is 
the simplest one and more complicated decay process 
will need higher derivative of the potentials. So, the 
quadratic potential model is stable if we consider only 
the tree level potential coupling. But the other two po- 
tentials need momentum cutoff of A < lO^^Afp, which is 
still acceptable Q. 

While, if we also consider the derivative coupling, 
which is potential independent, the situation will be the 
same as that in Ref. 14] and the optimistic estimation of 
the momentum cutoff for the models should be less than 
10^°Mp ~ lOOMeV. However, this is not necessarily im- 
possible for phantom if we assume the scalar field theory 
is only valid up to relatively low momenta or find some 
mechanism to suppress derivative couplings, leaving only 
the couplings of gravitons to the potential, which were 
consistent with a cutoff as high as the Planck scale p^. 
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The point is that if the observations eventually do sug- 
gest dark energy with an equation of state w < —1 and 
if it is modelled by phantom, then it must be stable with 
some mechanisms such as those described above. 



Using Ea. (|16|) . we find that the remaining time tr = 
tend — to before the end of the Universe is given ana- 
lytically by 



VI. TWO ANALYTICAL MODELS 



In the following, to expose the above discussion more 
specifically, we discuss two solvable Phantom models that 
correspond to quasi-de Sitter and big rip late time phase 
respectively. We consider only the case that Phantom 
become dominant. The potential for the model corre- 
sponding to quasi-de Sitter solution is in the periodic 
potential Ea. H12|l . For simplicity, we consider only the 
behavior of for (f> e [0o,</'o + '^^]- The quasi-de 

Sitter solutions are 



a = ao{sec[ViK{(/) — 



(14) 



-{In 



(288Uo^)^«^ Vl-[l-e)«]^ 



+2arctan[l- (^)6]3} 

V2yo[i-(-)^] 

a 

1 



1 - (ao/ay 

It is easily found that the scale factor a{t) tends to infinite 
as — !■ 00+ 2v13k ' Furthermore, as the Universe expands, 
the Phantom energy density tends to ^/2Vo and the 
equation of state approaches —1. 

The second potential for the model admitting big rip 
solution is 



F((/)) =yoexp[v^AtA((/.-(/.o)] 



(15) 



where Vq and A are positive constant. The solutions of 
the system are 



Oq exp 



2(2 + ^2)1/2 
\/3k(2Vo)i/2A2 ^' 



(16) 



2Vn 



W4, = 



2 + A^^aa 
-(1 + ^2) 



Pol — 



(-) 

a 



y/^K.A'^ 



(17) 



In Ea. l(T7|l . putting in k 



V^Q.a,Hq with VIa, = 0.73 



and H, 







li.lGyr, one finds tr = 28.01, 53.4 and 106.9 



Gyr. for w = —1.38, -1.20 and -1.10 respectively. 
VII. DISCUSSION 



We study the phase space of phantom model in an ar- 
bitrary potential and give the corresponding conditions 
for tracking attractor, big rip attractor and de Sitter at- 
tractor. The big rip attractor and de Sitter attractor 
are two late time attractors that will determine the type 
of the cosmic doomsday. For the de Sitter attractor, 
we have A = and therefore a divergent F according 
to its definition. However, in the equations of motion 
Eqs.Q, such singularity will not appear because of the 
term A2r = ^ gives regular value. On the other hand, 
the tracking attractor serves a similar function as it does 
in traditional quintessence models and could alleviate the 
fine tuning of the model. The observational constraint 
— 1.38 < w < 0.82 is often interpreted as an observa- 
tional support for the existence of Phantom in our Uni- 
verse. We show that there maybe two possible phantom 
cosmic doomsday, big rip and de Sitter, which are both 
attractor phases of the dynamical system and compatible 
with current observations. At least before the w < —1 
dark energy is excluded by observation completely, the 
phantom models will remain an interesting alternative. 
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